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Partial eruption of a filament with twisting nonuniform fields 

Yi Yunchun Jiang\ Jiayan Yang\ Yongyuan Xiang^, Yunfang Cai^ AND Weiwei Liu^ 

ABSTRACT 

The eruption of the filament with the kink fashion is often regarded as a 
signature of the kink instability. However, the kink instability threshold for the 
filament magnetic structure has been not widely understood. Using the Ha ob¬ 
servation from the New Vacuum Solar Telescope (NVST), we present a partial 
eruptive filament. In the eruption, a filament thread appeared to split from the 
middle portion of the filament and to break out in a kinklike fashion. During 
this period, the left filament material remained below, which erupted without 
the kinking motion later on. The coronal magnetic field lines associated with the 
filament are obtained from the nonlinear force-free field (NLFFF) extrapolations 
using the 12 minutes cadence vector magnetograms of the Helioseismic and Mag¬ 
netic Imager (HMI) on board the Solar Dynamic Observatory (SDO). We studied 
the extrapolated field lines passing through the magnetic dips that are in good 
agreement with the observed filament. The field lines are non-uniformly twisted 
and appear to be made up by two twisted fiux ropes winding about each other. 

One of them has higher twist than the other, and the highly twisted one has its 
dips aligned with the kinking eruptive thread at the beginning of its eruption. 
Before the eruption, moreover, the highly twisted flux rope was found to expand 
with the approximately constant field twist. In addition, the helicity flux maps 
deduced from the HMI magnetograms show that some helicity is injected into 
the overlying magnetic arcade, but no significant helicity is injected into the flux 
ropes. Accordingly, we suggest that the highly twisted flux rope became kink 
unstable when the instability threshold declined with the expansion of the flux 
rope. 
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Introduction 


The kink instability is often considered a candidate mechanism for the ernption of 
the hlament. When the kink instability sets in, the flux rope’s axis will experience a 
writhing motion due to the conservation of magnetic helicity in the highly conducting corona. 


) 

ance of the kink instability 

(Rust & LaBonte 

2005; 

Green et . 

2007: 

Liu 

2008; 

Bi et al. 

2012; 

Thompson et al. 

2012; 

Yang et al. 

I 2 OI 2 ; 

Jiang et al. 

2 OI 3 I: ^ 

^an et al. ! 

2014 

1). Diherent 


from the rotation that is guided by s urroundirig ma gnetic held in the higher corona (e.g., 


Cohen et ahl 1201 Ol: iKliem et ahl 120121: iBi et ahl 120131 1 , the writhing motion induced by the 
kink inst ability mainly acts in the lower corona, below a height comparable to the footpoint 
distance (IKliem et al.ll2012l) . 


The ideal stability of the kink mode is mainly controlled by the total twist. The kink 
instability occurs if the amount of magnetic twist exceeds a critical value. The hux tube 
twist is usually expressed as ^tw = where B^/Br is the ratio of the azimuthal 

and axial held components of the hux tube and //r is the length-to-width ratio of the tube. 
The threshold for the onset of instability is dependant on the detailed mag netic structure. 
For a force-free magnetic loop including the effect of photospheric line tying (IHood fc Priest 
198111. t he ins tability threshold is an amount of twist about 1.25 turns. The MHD simulation 
of iFanI (1200511 showed that a line-tied hux rope becomes kink instability when the held line 
twist reaches about 1.7 turns. Consistently, the critical average twist for a magnetic loop 


simulated bv lTbrbk et ahl (120041 1 is about 1.75 turns. Moreover, the authors concluded that 
the instability threshold rises with the rising aspect ratio of the loop. 


An inhomogeneous twisted hux rope was investigated by iBirn et al.l (120061 1 using mag- 
netohydrodynamic simulations. Their results show that a hux rope will be broken into two 
portions when the more strongly twisted portion becomes kink-unstable and i noves rapidly 
outward. The bifurcation of a hux rope during er uption was hrst discu ssed bv ICilbert et al. 


(120001 1 and was classed as a partial eruption by ICibson fc FanI (120061 1 , who demonstrated 
that part of the hux rope is expelled from t he corona when it rec onnects internally and with 
surrounding held so that it breaks in two (iGibson fc Fanll2006l ). Their simulation implies 
that the writhing motion induced by the kink instability is essential for forming a current 
sheet within the hux rope where it can break in two. The observ ations of the partial eruption 
(e.g., Tripathi et al. 2009 : Shen et al. 2012 : Tripathi et al.l[201^ exhibit that the splitting of 
the hux rope is often accor npanied by the writhi ng motion. The authors accounted for these 
events using the model of iGibson fc FanI (120061 1. but it seems that not all of these events 
show a clear signal of the internal magnetic reconnection. 


The twisted hux rope is often found above the magnetic neutral line in the nonlinear 
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force free magnetic field (NLFFF) extrapolated from vector magnetograms. Using various 
NLFFF extrapolation algorithms, several au thors found the flux rope systems in t heir studied 
active region with tw ists of about 1 turn ( Yan et ^ boOl : Canon et ah 2009 : Inoue et ah 


2011 


Guo et al. 


■Ting et al.ll2ni4 


20131. Some filam ents investigated ( Regnier fc Amari 20041: IGuo et ah 2010 : 


■Tiang et al.ll2014Jl were also described by the twisted flux rope with the value 


of the twist of approxim ate 1 turn, w’hich mean that these structures are stable against kink 
instability. Interestingly, iGuo et ahl (120131 1 reconstructed a flux rope, which has the twist of 
about 1.9 turn about two hours before it started to erupt. It is worth noting that the ratio 
of the radius to the length of the flux r ope they stiidied is less than 1/15, which is slightly 
lower than the ratio that was chosen bv lTbrdk et ahl (12004 1 in their parameter investigation. 
To our knowledge, however, few investigation has focused on the evolution of the field twist 
associated with a filament showing the kinking motion. 

Magnetic helicity quantifies the magnetic topological complexity and is a valid tool for 
measuring how much a magnetic flu x rope is twisted and writhed, or how much a magnetic 
arcade is sheared (see the review by iDemoulinI (120071 1 and references therein). Inferred from 
the evolution of the photospheric magnetic held, magnetic helicity transported across the 
photosphere is often used to diagnose the helicity that is instantaneously injected into the 
higher solar atmosphere. Some investigates showed that the local injection of helicity with 


oppo s ite signs played a r ol e in triggering the erup tions of the studied filaments (IRomano et al. 


201 ll : iDhara et al.ll20l4 l. iRomano et ahl ([20051) found the impulsive input of helicity at the 
beginning of the eruption of a kinked filament. Their results supported the idea that the 
amount of the magnetic helicity in the filament exceeds the limit for the kink instability 
primarily due to the transport of the helicity through the photosphere. 


In this article, we present the partial eruption of a filament. One portion of the filament 
shows the evident kinking motion in the eruption. We present observations of the event and 
comparisons with the topologies of the NLFFF extrapolated held lines and the photospheric 
helicity injection associated with the filament. This enables us to investigate the structure 
and the evolution of the magnetic field relating to the kinking eruption. 


2. Observation and data analysis 
2.1. NVST and SDO observations 


Our primary data, t o be used to present the eruption of the filament, are the Ho 
images obtain by NVST ( Liu fc Beckers 2001: iLiu et al. 2014 1 and the 304 A images from 


the Atmospheric Imagining Assembly (AIA; 


Lemen et al. 


20121 1 on board SDO. The NVST 
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Fig. 1.— (a) NVST Ha image, (b) A line-of -sight (LoS) magnetogram in the CCD 

coordination, (c) The radial component of the vector data in the CEA coordination. The 
FOV of panel (a) and (b) is the full view of the NVST and the SHARP region, respectively. 
On panel (b), the dash box colored blue indicates the position and FOV of panel (a) and 
the red box indicates that of panel(c) The solid blue boxes on the panel (b) and (c) show 
the same position. 


obtains the solar images in the three channels, i.e., Ha, TiO, and G band, with a field of view 
(FOV) of 180" X 180", a cadence of 12s, and a pixel size of 0".168. The observational data 
acquired at NVST are now available (http://fso.ynao.ac.cn). The NVST Ha adopted here 
were obtained from 04:00 UT to 07:00 UT on 2014 Nov 4, and the whole FOV of NVST/Ha 
image is displayed in Figure 1(a). All of the NVST images are aligned to the same FOV 


based on a high accuracy solar image registration procedure (IFeng et ahl 120121 : lYang et ah 


20141 ). The AIA takes full-disk EUV images with a pixel size of 0".6 and 12s cadence. 


The HMI vector field (ITurmon et al.ll2010ll is computed using the Very Fast Inversion 


of the Stokes Vector IVFISV: iBorrero et al.ll201ll ) code and the remaining 180° azimiith arn - 
biguity is resolved with the Minimum Energy (MEO) code (lMetcalllll994l : iLeka et ahl 120091 ). 
Now, HMI provides continuous coverage of the vector field on the so-called HMI Active Re¬ 
gion Patches (HARPs) region. Figure 1(b) present a line-of-sight (LoS) magnetogram in the 
CCD coordination, with the same FOV as the HARP region surrounding the targeted fila¬ 
ment. The HARP vector field data has been remapped to a Lambert Cylindrical Equal-Area 
(CEA) projection and then transformed into a standard heliographic spherical coordinates. 
Figure 1(c) displays the CEA map-projected field in the radial direction normal to the solar 
surface. 
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2.2. Coronal magnetic extrapolation 


The “NLFFF” package available in SSW is develope d by Jim McTiernan to perform a 


Wheatland et 


NLFFF extrapolation using the optimization method of 
one o f the best-performing NLFF field modeling methods fISchrijver et ah 


al. 


nn^ : 


( 20001. which is 


Metcalf et al. 


20081 1. The procedu re is fulfilled in the Cartesian and spherical coordinate and includes the 
weighting function (IWiegelmannI 1200411 and a preprocessing proce dure to drive the observe d 
data towards suitable boundary conditions for an extrapolation (jWiegelmann et al.l 1200611 . 
In this study, the coronal field on the spherical geometry is extrapolated by means of the 
“NLFFF” package as follows. To ensure the vector dat a for use by NLFFF modeling efforts 
spans the area as large as possible flDe Rosa et al.ll20091 1 , firstly, the field covering a relatively 
large area of 21° x 21° (indicated by the dashed box in Figure 1(b)) is calculated within a 
computational domain of 200 x 200 x 240 grids. The CEA vector field is used as the bound¬ 
ary condition and the resul ts of the potential-field source surface (PFSS dSchatten et ahlllOfiOl: 
Schrijver fc De Rosall2003l ) model as the initial condition for this process. Secondly, the so¬ 
lution field is applied as the initial conditions for the further extrapolation within a domain 
covering an area of 16° x 9° that includes the region of interest (indicated by the solid box 
in Figure 1(b)). The domain is resolved by 534 x 301 x 477 grids to ensure that the size 
of the grid is same as the pixel size of the HMI data. Fo r evaluating the perform ance of the 
NLFFF extrapolation, the (dj) metric and (|/i|) metric flWheatland et ahl 120001 1 range from 
16° to 22° and from 6.7 x 10“^ to 7.5 x 10“^, respectively. The value is enough close to zero 
to ensure that the field is close to the force-free and divergence-free state. 


2.3. Twist and writhe of the magnetic field line 


A modeled field line may indicate a thin flux tube when it is integrated with its sur¬ 
rounding almost parallel field. The helicity of a tube can be decomposed into contributions 
from twist and writhe, i.e. 

H = (1) 

where $ is the magnetic fl ux of the tube, the v alues of T^, and Wr are deduced from the 
geometry of the field lines flBerger fc Priori 1200611 . such as 




( 2 ) 


and 


w,= ^<f T(s) X T(s') ■ 


( 3 ) 
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in which the the tangent(T) and normal(N) unit vector satisfy the Frenet-Serret formulas. 

For a uniform flux rope with the twist ^tw = (see Figure 2), is ^{IB^/r^B^ + B‘^) 

according to the Equation (2). As B^/B^ is often smaller than 1 in the solar corona, the 
value of Tw is close to that of defined in the theoretical calculation. In this study, 

Tw is used as the proxy for the twists of the modeled field lines. 


B 



B 



Fig. 2.— Sketch of two field lines, which were generated using the form: x(t) = B^sin{t)-, 
y(t) — rB^cos(t)\ z{t) = Brt with t G [0,l/br], where l/r — 20, B^/Br = 3/4, so that the 
twist = IB^/rBr = 2.47r, and IB^/r^B"^ + B^ = l.Ovr. According to the Equation (2), 

Tyj for the field lines on the left (resp. right) is 1.9 (resp. -1.9). The unit vectors T, N, 
and B, often called the tangent, normal, and binormal unit vectors, satisfy the Frenet-Serret 
formulas. 


2.4. Magnetic helicity injection 


The rate of helicity transport across th e photosphere is often computed via the proxy 
0$ fjPariat et al.ll2005l : iDalmasse et al.ll20141 1. 0$ is a function associated with the magnetic 
flux transport velocity, which is determined by local optical flow technol ogy, such as l ocal 
correlation tracking (LCT)or differential affine velocity estimator (DAVE; ISchnckI 120951) . In 
this study, DAVE is implemented with the window size of 20 x 20 pixels and the time 
interval of 720s. 
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It is worth noting that the gauge invariance of the helicity flux depends on integrating 
over the surfaces or at least a big patch. To b etter estimate the spatial distribution of 
the helicity flux, accordingly, IPariat et al.l (l2005[ i considered the helicity injected into an 
elementary flux tube, which can be written as the sum of the Gg at the both photospheric 
footpoints of the corona connection. The authors introduced a connectivity-based helicity 
flux density proxy, G^, such that 


G(f){x±) — -I- Gs^Xzp) 


Bn{x±) 


Bn{Xzf) 


f± 


(4) 


where x+ and X- are the locations of the opposite polarities with coronal linkage, and f 
ranges from 0 to 1. Then proxy G^ has a meaning for the helicity flux for each individual 
flux tube. Obviously, the integral of G^ has the same value of that of Gq over the region 
large enough, but G^ masks some false opposite sign signals in Gq, which is canceled out in 
the integral of Gg over the region including all magnetic polarities with coronal linkage. 


3. Results 


The close views of the filament on the Ha and 304 A image are shown in Figure 3. 
The filament is located in the NOAA active region AR11884 at a position (Sll, W24) on 
4 October 2013. The length of the filame nt is approximately 30 Mm, slightly longer than 
that of the named mini-fllament defined by Wang et al. f 2000li . At 05:28 UT, a thin filament 
thread starts to split from the middle portion of the filament (as indicated by the arrows in 
Figure 3(c-d)). The separated thread was found to rotate clockwise. Observed from the Ha 
image(Figure 3(c)), the rotation angle about the undisturbed filament axis is approximately 
30° at 05:36 UT (also see the accompanying animations of Figure 3). At 05:38 UT, the 
eruptive filament thread has been disappeared from the Ha image (Figure 3(e)) and it 
evolved into a brightened one in the AIA 304 A image (Figure 3(f)). The left portion of the 
filament remained undisturbed until after the previous eruptive thread has disappeared in 
the AIA 304 A image. Moreover, no evident rotating motion is identified in the subsequent 
eruption of the left portion (Figure 3(g-h)), which started about 05:45 UT. Ultimately, a 
flare of X-ray class C3.2 took place with start, peak, and end times around 05:36, 05:44, and 
05:52UT, respectively. 


The modeled magnetic dips exist if the held lines are locally horizontal and curved 
upwardly, and they are widely accepted to support the filament plasma. In Figure 4(a) and 
4(b), the green curves overlying the Ha images show the magnetic dips calculated from the 
extrapolated NLFFF field at 04:12 UT and 05:12 UT, respectively. At these two moments, 
the locations of the dips seem to be very well aligned with the filament over its length. 
















Fig. 3.— Left: NVST Ha image. Right: SDO AIA 304 A images. 

The held lines traced from all of the dips are displayed in the (p — R and (p — 9 planes 
(Figure 5). After checking the character of the each held line, we assign the held lines into 
two groups: the held lines are colored red if their endpoints (as showed in Figure 4(c-d)) that 
anchor in negative polarity have the held strength higher than 1200 G, otherwise, the held 
lines are colored cyan. After such a grouping, the held lines seem to be made up by two hux 
ropes winding with each other(see the First two rows of Figure 5). Red colored held lines 
show a double peaks structure with their two peaks being located above the most of cyan 
colored held lines, while the majority of cyan held lines have their peaks lying between the 
two peaks of red held lines. Moreover, the dips of the red colored held lines match well the 
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04:12 UT 



28 29 30 31 

Lontitude (Degree) 


05:12 UT 


28 29 30 31 

Lontitude (Degree) 



Fig. 4.— (a-b): The NVST/Ha images overlaid by all of the extrapolated magnetic dips, 
(c-d): The radial component of the HMI vector field overlaid by the endpoints of the fields 
traced from the location of the dips. The endpoints are colored red if their corresponding 
fields anchor the negative region having the field strength higher than 1200 G, otherwise, 
they are colored cyan, (e-f): The NVST/H« images overlaid by the extrapolated magnetic 
dips, which are associated with the fields that anchor the negative region having the field 
strength higher than 1200 G. All of these images are remaped to the CEA coordination. 
(Animations of this figure are available in the online journal.) 
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04:12 UT 


05:12 UT 




28 29 30 31 

Lontitude (Degree) 


28 29 30 31 

Lontitude (Degree) 



28 29 30 31 

Lontitude (Degree) 


28 29 30 31 
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Fig. 5.— The field lines traced from the magnetic dips as showed in Figure 4(a-b) are 
displayed in the (j) — R and (p — O plane. The fields are colored red if they anchor the negative 
region having the field strength higher than 1200 G, otherwise, they are colored cyan. 

middle segment of the observed filament (Figure 4(e)). In figure 4(f), the eruptive filament 
thread appears to be lying across the middle part of the filament, where approximately the 
dips are located. 
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After obtaining the valne of the twist and the writhe for each field line according to 
the equation (2), we found that the extrapolated field fines are non-uniformly twisted. The 
values of the twists range from 0.8 to 2.0 turn and average about 1.2 turn. Furthermore, 
the average twist for each field fine colored red is higher than the average for the field line 
colored cyan. The averaged value for the red (resp. cyan) is 1.55 (resp. 1.25) at 04:12 UT 
and 1.51 (resp. 1.1) at 05:12 UT. The value of the writhe is smaller than that of the twist. 
The average writhe of the field fines colored by red has the value of -0.04 at 04:12 UT and 
-0.02 at 05:12 UT. 

Possessing higher twist, the flux rope colored red is extracted to be further concerned 
on the last two rows of Figure 5. The flux rope has the height of about 7 Mm referred to 
the solar surface, and has the length of about 30 Mm. Comparing the two panels, we can 
find that the flux rope evolves to be fatter. The averaged radius of the flux rope is about 2 
Mm (resp. 4 Mm) at 04:12 (resp. 05:12) UT. Hence, the ratio between the radius and the 
length of each flux ropes is about 1/15 and 1/8 at these two moments, respectively. 



29 30 31 29 30 31 29 30 31 

Lontitude (Degree) Lontitude (Degree) Lontitude (Degree) 


Fig. 6.— (a) The radial component of the HMI vector field at 05:12 UT. (b) The time- 
averaged Gg map. (c) The time-averaged map. On panel (a), the selected field lines 
connecting the N2 are colored blue, and the endpoints of the field lines as the Figure 4(d) 
showed are overlaid. On panel (b-c), the color shows the strength of the helicity flux with 
white being 0, red positive, and blue negative. The radial component of the HMI vector field 
is contoured on these panels. The contour levels are 40G for positive polarities (solid) and 
-40G for negative polarities (dash). 

Figure 6(b) displays the time-averaged helicity flux map computed with the proxy Gg 
by a series of HMI magnetograms at 12 minute cadence from 04:12 UT to 05:36 UT. Figure 
6 (a) exhibits one of the magnetograms, on which the two negative polarities is referred as N1 
and N2, respectively. From the NLFFF, we note that the N1 is where the eastern endpoint 
of the filament anchors and the N2 connects the field overlying the filament, while no other 
negative patch flux is associated with the filament. And yet the part of the positive polarity 
around the western footpoint of the filament connects the large-scale field that anchored far 
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away from the filament. Accordingly, we choose f=l to compute the G^{x-) (Figure 6(c)), 
which assumes the helicity is injected from negative polarity. In the G^{x_) map, therefore, 
helicity flux from the N1 and N2 are roughly equivalent to the total helicity injection into 
the filament. 

As showed in Figure 6(b-c), both the Ge and G^ maps display mixed signals in Nl. 
However, the G^ map presents the more positive flux on the right side in Nl, and then the 
total helicity flux in Nl of G^ is close to 0. It is consistent with the result deduced from the 
extrapolated field lines, which show that no increase of the magnetic helicity was discovered 
in the magnetic structure associated with the filament. In contrast, the distributions of Ge 
and Gcf) are similar in N2, which is completely filled with positive helicity flux. 


Conclusions and discussion 


By means of observations from NVST and SDO, we present a partial filament eruption. 
The eruption starts with a thin filament thread being separated from the filament, and 
the thin thread breaks out in a kinklike fashion. We regarded the field lines with their 
dips aligning the observed filament as the magnetic structure associated with the filament, 
and investigated the twists and writhes of the field lines based on their geometries. The 
result demonstrates that the magnetic system consists of two flux ropes, and one of them 
is estimated to be more highly twisted than the other. Moreover, the highly twisted flux 
rope has its dips aligned well with the thin filament thread that shows kinking motion in 
eruption. Therefore, the splitting of the eruption may result f rom that the mo re strongly 
twisted portion of the magnetic system becomes kink-unstable (IBirn et ahl 120061 1. 


For each field line constituting the highly twisted flux rope, the average twist is found 
to decrease slightly with a small increase of the average writhe as time passes. Hence, some 
twist may have been quasi-statically converted into the writhe duri ng this per iod. A similar 
equilibrium configuration was obtained in the MHD simulations of iFanI (120051 ). who showed 
a mildly kinked equilibrium with some finite writhing of the flux rope axis. Moreover, the 
helicity for the highly twisted flux rope, which is equivalent to the sum of the values of the 
twist and the writhe according to the Equation (1), amounts to 1.51$^ and 1.49$^ at 04:12 
UT and 05:12 UT, respectively. It indicates that no increase of the amount of the magnetic 
helicity was detected in the flux rope before its eruption. In this case, the question is the way 
in which the flux rope develops to be nonlinear kink instability as implied by the observed 
strong rotation of the filament thread. 


As mentioned in the Section 2, the value of is slightly smaller than that of the twist 
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^tw = IB^/rBj. for the same flux rope. It implies that the twist ^tw for our mode led flux rope 


should be slightly greater than 1.5 turns that is deduced from the proxy Tyj. iTorok et al. 


(120041 1 estimated ^tw = 1.75 turn as the threshold of the kink-instability of a magnetic loop 
with their defined loop aspect ratio of about 5, which amounts to a ratio of 11.2 between the 
width and the length of the loop. In our extrapolation, the width-length ratio of the modeled 
flux rope incr ease from 1/ 1 5 to 1 /8 before the eruption, i.e. from less than to greater than 
that given by Tor ok et al. ( 2004ll . As the instability threshold rises with rising aspect ratio 
(ITorok et ahl 120041 1 . it is possible that the threshold for the kink instability is greater than 
1.75 turn at 04:12 UT and is less than 1.75 turn at 05:12 UT in our extrapolation. Hence, the 
field structure may become kink unstable when it expands with the approximate constant 
of the twist. 


On the other hand, the helicity flux from the photosphere is obtained using both the 
proxies Gq and G^. Since the proxy G(f, involves field connectivity of each elementary flux 
tube, the proxy G^ is better than the proxy Gq to present the distribution of photospheric 
helicity flux density per flux tube. However, direct interpretation of the distribution of the 
helicity flux should be made with caution, because the determination of the helicity flux 
requires the gauge selection for the definition of the vector potential A(V x A = B) and 
both the helicity proxies Ge and G^ are calculated with the classical gauge of V x Ap — 0 
and Ap ■ n = 0. By contrast, there is no gauge involved in T^, and H7 according to the 
Equation (2) and (3). In our study, the results from both methods show that no significant 
magnetic helicity is injected into the magnetic flux rope passing through the magnetic dips. 


Moreover, both proxies Go and G^j^ exhibit that some positive helicity is injected from 
a patch flux (referred as N2) which connects the magnetic arcade overlying the studied 
filament. The local injection of the helicity may play a role in twisti ng the overly i ng arc ades. 
If the overlying field is weakened by the twisting, as suggested bv ITorok et ahl (j2013[ (. the 
underlying flux rope would undergo a quasi-static expansion. Therefore, the helicity injection 
into the overlying arcade may have a role in expanding the flux rope as showed in the 
extrapolation. 
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